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We adopt methods from statistical field theory to stochastic inflation. For the example of a free 
test field in de Sitter and power-law inflation, the power spectrum of long-wavelength fluctuations 
is computed. We study its dependence on the shape of the filter that separates long from short 
wavelength modes. While for filters with infinite support the phenomenon of dimensional reductions 
is found on large super-horizon scales, filters with compact support return a scale-invariant power 
spectrum in the infra-red. Features of the power spectrum, induced by the filter, decay within a few 
e-foldings. Thus the late-time power spectrum is independent of the filter details. 
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I. INTRODUCTION 

Despite of its successes as a building block of our cur- 
rent picture of the Universe, a full understanding of the 
inflationary dynamics on super-horizon scales is still lack- 
ing. In his pioneering work yQ, Starobinsky introduced 
the concept of stochastic inflation to provide a frame- 
work to study the evolution of quantum fields in an 
inflationary universe [2J [3J HJ [5J [7J. The key idea, 
which acquired considerable interest over the last years 
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lies in splitting the quantum fields into long- and short- 
wavelength modes, and viewing the former as classical 
objects evolving stochastically in an environment pro- 
vided by quantum fluctuations of shorter wavelengths. 
Hence, it constitutes an example of how fundamental 
properties of quantum fields can be modeled using meth- 
ods of statistical mechanics. 

Given the de Sitter horizon, \, as a natural length 
scale of the problem, one then focusses on the "relevant" 
degrees of freedom (the long-wavelength modes) and re- 
gards the short-wavelength modes as "irrelevant" ones, 
where "short" and "long" are subject to \- 

The most simple setup provides a fixed cosmological 
background, in which the dynamics of a scalar test field 
tp is analyzed. If tp is free, massive and minimally cou- 
pled, one obtains after splitting into long and short wave- 
lengths, ip = <y9 L + ip s , an effective equation of motion of 
generalized Langevin-type, 



p + fj, 2 )<p h (t,x) =h(t,x) 



(1) Sec 



In equation Q, <^ L is viewed as a classical entity, evolv- 
ing stochastically in the presence of a (quantum) random 
force h, which is Gaussian distributed with zero mean. 

Self-interactions cause deviations from the simple 
Langevin-type form, manifested in higher powers of ip L 
on both sides of Eq. ([!]) with coefficients being of stochas- 
tic origin. The methods presented in this work allow to 



analyze this most general test-field case, which, to our 
knowledge, has not been addressed so far. Early studies 
focussed on homogeneous fields, thus restricting atten- 
tion to the time evolution of y L . Recently, we presented 
a method to calculate arbitrary 2-point functions for gen- 
eral stochastic potentials [22] . 

This work extends our previous results [22] and is de- 
voted to the study of the scaling behavior and time evolu- 
tion of the power spectrum of <p L in a fixed background. 
This is done by means of replica field theory, which is 
well known in statistical physics [23 . Replica field the- 
ory allows us to compute the non-coincident two-point 
function of the long- wavelength modes for the most gen- 
eral test-field case. In order to distinguish short from long 
wavelength, a filter is introduced. The most common fil- 
ter in stochastic inflation is a (sharp) step function. We 
study in detail the dependence of the power spectrum 
on the shape and parameters of smooth filter functions 
and consider the sharp step as a limiting case. For filter 
functions with infinite support, a variant of the so-called 
dimensional reduction [21 [251 [251 [27J [251 is found, 
which results in a strong deviation from scale invariancc 
in the infra-red, signalling a breakdown of the test-field 
assumption [22] . However, we show in this work that 
scale invariance is preserved on all scales at late times 
for filters with compact support. 

Our work is organized as follows: After calculating the 
free, noiseless propagator for exponential as well as for 
power-law inflation (Sec. [H]), we add quantum noise in 
Sec. |TTX Its effect on the power spectrum is discussed in 
Special focus is put on the dependence on the 
filter functions and a non-linearity parameter g NL is in- 
troduced to quantify the modification of the power spec- 
trum. The methods needed, are described in appendices 
[A] and [5] In particular, the replica trick is introduced and 
a variational technique is presented. We conclude with a 
summary. 



II. FREE POWER SPECTRUM 
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This section is devoted to review the quantization of 
a free, minimally coupled, iV-component, real test field 
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ip with mass /z. We solve the classical mode equations, 
construct the propagator, the power spectrum and give 
expressions for the spectral index. 

Let us concentrate on a spatially-flat, isotropic and ho- 
mogeneous universe in four-dimensional space-time. For 
its scale factor we assume cither exponential inflation, 
a{t) := e Ht , with Hubble rate H, or power-law inflation, 
a(t) := (t/ti) p with p > 1 and the reference time t\ de- 
fined by a\ '■= a(ti) = 1. For convenience use H = c = 1 
and set either H — 1 or ii = 1, respectively. The mode 
function u(t, k) is defined via the decomposition of the 
field components (i = 1 , . . . , N) 



<Pi(t, k) = &i(k)u(t, k) + H.c. 



(2) 



with the modulus of the comoving momentum k := |fe|. 
The annihilation and creation operators obey the com- 
mutation relations 

[& j (fc),&t(p)] = (27r) 3 d 3 (k-p)5 lJ1 (3a) 

[&*(*;),%(?>)] =0. (3b) 

The rescaled mode functions v(t, k) := a(r) u(r, fc) ful- 
fil the mode equation 



,2 , 2 2 

k + (i a — 



v = 0, 



(4) 



with primes denoting derivatives w.r.t. conformal time 
r, defined by dr := a(i) -1 di. Solutions to Q are fixed 
by requiring that for very short wavelengths the effect of 
space-time curvature and mass becomes irrelevant, and 
thus a plain-wave solution should be obtained, i.e., 



lim v(t, k) = 



~— i k 



12k 



(5) 



The factor l/v2 k is fixed by the canonical commutation 
relations of ip and its conjugate momentum. 

Exponential inflation implies a(r) = |r| 1 and for 

_ p 

power-law models one finds a(r) oc \t\ p_1 , where both 
cases match in the limit p —> 
be rewritten in the form 



Then, equation Q can 
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v = 0, 



(6) 



with 



/' = i/ - -/r 



(exponential), (7a) 



and with /i = 



» 1 

z/ = 1 — (power-law). (7b) 

p—1 2 

A general solution to |6]) , fulfilling ^ , is given in terms 
of Bessel functions: 



,1/2 



m(t, fc) 



2 a(r) 



(J^fclrD+iY^fcM)). (8) 



On large scales (for k — > 0) and for positive the leading 
term of ([8]) is 



u(t, k) ~ — 



i2 *-i i T .i/a- 



/7ra(r) 



(9a) 



while for negative v one has 



-v-\ 



1 1/2+1/ 



u(r, fc) 



(i COt(7T i/) 



r(i/ + i)o(T) 

The propagator is defined as 

G Q (tXk,k') := l(fi|^(<,fc).^',fc')|"), 



r. (9b) 



(10) 



where the vacuum is defined by a(fc)|f2) = and 
a subscript "0" indicates a quantity that is calculated 
in the absence of any classical noise. The factor 1/N is 
convenience. 

An object of central interest in cosmology is the di- 
mensionless power spectrum r P v (k). Its relation to some 
field propagator G(k), with assumed infra-red behavior 
G(fc) ~ k~ v , is given by 



Vip(k) := k 3 G(k) ~ k n *- 



(11) 



with the spectral index n v . It is connected to the so-called 
critical exponent rj via n v = 4 — r\. 

For n = the power spectrum of the free, noiseless 
theory is scale-invariant, i.e., n v = 1 for ryo = 3. Non- 
zero mass leads to 



3a/1 



1 + T u 2 + 0(n 4 



9 r 3 r 
and the power-law case with /i = yields 



p- 



1 



1 



2 
P 



O 



1 



(12a) 



(12b) 



Note, that n v (/i ^ 0) > 1 for exponential inflation, while 
riip(p < oo) < 1 in the massless power-law case, but in the 
limit [i — > 0, or, p — > oo, respectively, one recovers a scale- 
invariance spectrum. Also note, that the results (12a I 



and (12b I do not include any metric perturbation, which 
would generally cause negative deviations from n v — I 
(see, e.g., [5U]). 



III. ADDING NOISE 

Here, we describe the dynamics of the long- wave length 
modes, <^ L , of some quantum field ip and how they 
are influenced by its short- wavelength counterparts. In 
this sense, the notion of stochastic inflation refers to a 
stochastic modelling of the former (treated as classical) 
evolving in a (stochastic) environment provided by the 
latter. Now, how does one obtain the aforementioned 
Langevin-type equation? 
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To demonstrate this idea, let us follow Kandrup [6 
[39j and consider (in four dimensions and for N = 1) the 
equation of motion 



(□ + n 2 )(fi(t,x) = 0. 



(13) 



Then one splits the field tp into a long and short- 
wavelengths part, (p = </? L + ip s , like 



<p s (t,x) = J d 3 k e(kr(t) - e) : 



a(fe) u(t, k) e 



— i fe-; 



a\k)u*(t,k)e 



k-x 



(14) 



is the usual step function that sharply cuts out the long 
wavelengths. Instead, one could also write an arbitrary, 
smooth filter function (c.f. the discussion in section IV C I. 



The parameter e, which basically says where to split 
into short and long wavelengths, should be chosen much 
smaller than one in order to cut far beyond the Hubble 
radius. On the other hand e is bounded from below, since 
otherwise, the change in the background metric has to be 
taken into account. This leads to e ^ exp{— H 2 /H}, i.e., 
|ln{e}|»0 
Inserting ( |14j ) in ( |l3| yields 

(15) 



(D + fi 2 )ip L (t,x) =h(t,x), 
with h being given by 



h(t,x) := / d 3 k5(kT(t) -e)x 



a(fc) w(t, k, x) e 



-i k-x 



aJ(fe) w*(t, k, x) e^ 



i k-x 



(16) 



Here, w(t, k, x) is obtained after the application of the 
covariant Laplacian (see [B] for details). The point now 
is, that <yS L is viewed as a classical entity (instead of as a 
quantum one) which evolves in the presence of the "ran- 
dom" force h that is itself treated as a genuine quantum 
object. This turns equation (151 into a classical (gener- 
alized) Langevin equation with h being a Gaussian dis- 
tributed random variable, with the moments (quantum 
averages) (h(t,x)) = and (h(i, x) h(t', x 1 )) ^ 0. 

For exponential inflation the noise two-point function 
is approximately given by [6] 



(h{t,x)h(t',x')) 



47T 



( e \x — x'\ \ 

9 sin 1 

r(t) 



However, this result heavily relies on the particular choice 
of the filter function. The Markov property (manifest 
through the time delta function) is a consequence of the 
step function and does in general not occur for arbitrary 
smooth window functions. A more detailed discussion on 
these functions is presented in section [IV C[ 

The Lagrangian £(y L ) for the long-modes <p L corre- 
sponding to equation (15 1 has, despite of its free part 



Co — associated with the l.h.s. of ( 15 1 — two additional 



pieces. The first provides the interaction with the noise 
and a second is quadratic in h, which is ascribed to 
a stochastic probability distribution p[h] [c.f. equation 
(21 1]. The Lagrangian (including only y> L ) is 



£(<£ L ) = £o(<^l) +h<yS L 



(18) 



to which, in statistical field theory, is referred to as the 
random-field (rf) case. As has been shown recently, this 
effectively re-sums the leading-log contribution of the full 
quantum theory [31]. Please note that equation (18 1 has 
been derived for free fields only. 

As soon as interactions are taken into account it can 
(and in general will) be only an approximation. Since 
one is often interested in theories containing interaction 
potentials, and in addition in situations where quantum 
fluctuations may become large, a generic effective inter- 
acting theory will generate some stochastic potential V, 



(19) 



meaning that the Taylor coefficients of V(<^ L ) are ran- 
dom variables (see below). Here and in the following, 
we generalize to iV-component fields, with correlations 
Gij (t, k) = dij G(i, k). To simplify notation, we will drop 
the subscript 'L' in the following. 

In the random-field case the potential V corresponds 
to V SF (0(x)) = — J2j fj{ x )i whereas the case 
of so-called random anisotropy (ra) is described by 
Vra((^») = - Y^jk^3k(x) (Pj(x) <p k (x). In general one 
may has 

oo N 

(20) 



n— 1 ii...i n —l 



In all cases {h} is assumed to be a set of Gaussian- 
distributed random variables which' distribution may be 
written as 



p[{h}] cx expi - i^h A A A ,.vhv - ^Mxh A 1 

^ A, A' A ' 



(21) 



with mean fi\. To simplify notation further, we use a 
general index containing position and all other compo- 
nents, i.e., we write h\ instead of h^...^ (x). 

Distribution (21 1 is an approximation for general inter- 



actions and justified for a limited class of theories only. 
However, for many systems it is at least a reasonable 
starting point, providing the possibility to obtain exact 
non-perturbative results. In this work we focus on the 
random- field noise, i.e., free, massive scalar test fields in 
various backgrounds. 

For some quantity O depending on {h}, the stochastic 
average, which shall be denoted by a bar, is then calcu- 
lated as 



0[{h}] := 2?[{h}]p[{h}]0[{h}]. 



(22) 
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Note, that linearizing the equation of motion in the quan- 
tum modes, corresponds exactly to a Gaussian distribu- 
tion of mean zero. Therefore, (211 is not an approxima- 
tion in this case. 

The methods needed for the treatment of these 
stochastic systems are discussed in appendix [A] and have 
been first introduced in the cosmological context in our 
recent work [2~21. 



IV. EFFECTIVE POWER SPECTRUM 

This section contains the study of the infra-red be- 
havior of the physical propagator G and therefore of the 
power spectrum V, with special focus on the role of the 
filter. 



A. Long-Range Correlation 



The full two-point function for the long- wavelen gth 

modes, G := jt(<Pl • <Pl,), is calculated in appendix B 
using a Feynman-Jensen-type variational calculation to- 
gether with replica field theory. The result is 



Equation (24) constitues one variant of the pheno- 



G(t, k) = G (t, k) + o(t, k) G (t, kf 



(23) 



where the so-called replica structure a is determined 
from the variational equations (B14al and (B14cl, re- 
spectively. 



To analyze the physical consequence of equation ( 23 ) 
on the power spectrum, let us now assume for the 
two-point noise correlation hi(x)hj(y) = <j>{ x — y) ~ 
\x — y\~ +p with p < 3 (c.f. appendix B2|. This is 
the case of so-called long-range correlation, introduced 
in [22] , and describes properly the infra-red limit of the 
physical model discussed below. In momentum space, 
the above choice implies 4>{k) ~ k~ p and hence a(t, k) = 
p(t)k~ p by virtue of (B14al. For p > —r)o, the infra- 



red behavior of the power spectrum deviates from the 
noiseless result, and we find r\ = 2?/o + P- This result is 
consistent with previous studies in flat space for a prop- 
agator with r; = 2 p?2"] . 

What does this mean for the spectral index? Since 
n v = 4 — r), we find the result 



n v = 4 - 2??o - p 



(24) 



if p > — T)q. Thus, one would find a dramatic change 
of the super-horizon power spectrum as compared to 
the case without noise. Concretely, for exponential in- 
flation this implies a modification of the spectral index 
on large super-horizon scales if the spatial noise corre- 
lator decreases at most like \x\ 6 if p = 0, while for fi- 
nite mass this exponent changes to —6 + 2/3 p 2 + 0(p i ). 
In the massless power-law model the power is given by 
-6-2p- 1 +C(p- 2 ). 



menon of dimensional reduction [26\ 128] . which can rig- 
orously be proven to all orders in perturbation theory for 
p = and for arbitrary non-random potentials (see es- 
pecially |28j for a supersymmetric version of the proof). 
Because this effect originates from the second piece of G 
in (23 1 it will be referred to as the dimensional reduction 
part. 

Please note, that those changes only concern the power 
spectrum of the smoothed (classical) long-wavelength 
modes, which are influenced by their short (quantum) 
counterparts. It does not mean that the full quantum 
two-point function obeys dimensional reduction. 

We should underline that the above statements on di- 
mensional reduction do depend on the concrete choice of 
the filter function. Their influence on the power spec- 
trum is discussed in section IIV CI 

A natural question is to ask on which scales does the 
effect of dimensional reduction show up. Let us therefore 
define the transition scale fc* at which the two terms on 



the r.h.s. of equation (23 1 balance each other via 



G (t, fe.) = a{t,K)G {t,Kf 



(25) 



It separates two regions such that for k 3> the behavior 
is noiseless and for k <C fc*. dimensional reduction holds. 



B. Stochastic Inflation 

Let us now return to our physical model of stochastic 
inflation. The split of the field <p into a long- and short- 
wavelength part, = (p^ + ips, together with the free field 
equation, (□ + p 2 )0 = 0, implies for the infra-red part of 
the propagator 



G(t,k) 



o(i) 



— e I u(t, k) 



\u(t,k)\ , 
(26) 



where is the (spatially Fourier transformed) covariant 
Laplacian, u(t, k) is the mode function from equati on (p i. 



and W K is a smooth high-pass filter (c.f. section IV C I, 
cutting out the low frequencies below e. The parameter 
k controls the width of the cut. In the limit k — > 0, W K 
approaches a step function. Here we choose 



W K ( • ) = — arctan 

7T \ K 



(27) 



and take < e < 1 in order to separate at wavelength 
well below the Hubble rate H(— 1), and n <C e to have a 
narrow transition region between quantum and classical 
modes. We do not impose any restr iction on p except 
that we demand the radicant in ( 12al to be positive, i.e., 
2 < 9/4. 



Please not that the filter function ( 27 1 does not have a 



compact support, meaning that also modes from the far 
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infra-red influence the quantum noise. A further discus- 
sion on filter functions, in particular of such with compact 
support, is presented in the subsequent section. 

Using ( 27 I the model given in ( 26 1 is of long-range-type 



(c.f. appendix LAI) and implies 



p = 3dl 



9 



2 = 1--M' 



+ 0(n 4 ) (28a) 



in the exponential case, and 



,P + 1 



1-p 



O 



(29b) 



for exponential inflation, and 
p+1 



P 



p-1 



i + 2 - + o fl 



P 



(28b) 



for massless power-law models. Thus for k <C fc* we 
obtain with ( 24 I 



6-9a/1 



-3 + 2^ 2 + (29a) 



for power-law inflation with /i = 0. We see that scale 
invariance of the effective power spectrum is destroyed 
even in the massless, exponential inflationary scenario 
in the infra-red limit. Also the power-law case changes 
drastically. 



For scales k~^> the noiseless spectral index (|12a 
recovered 

fc* , defined in equation ( 25 I 
cally: 



is 



The late-time behavior of the transition scale 
can be calculated analyti- 



8-2V9-4 M 2 

fc* = (e~*) n 



2\/9-4n 2 -2 



(e 2 + K 2 ) 



(30a) 



for exponential inflation and 



t p+i 



i i_ 

2 \ 2 2p 



7T 2 (e 2 + K 2 ) 



(30b) 



for massless power-law inflation. In the zero-mass limit 
(30a I yields the asymptotic form 



-t/2_ 



M 2 



(31) 



Thus, for e ^ 0, fc* goes to zero in the (step- function) 
limit k — » 0, i.e., dimensional reduction is absent — a 
statement that holds for exponential and power-law infla- 
tion as well. This is a general feature of the free (Gaus- 
sian) theory where is no mixing of the short quantum 
modes with the long classical ones as a sharp cut-off is 
introduced. In a slightly different setup, with a filter 
function having only one parameter, i.e., e = k, this has 
already been noted in reference |33j . 



Let us turn to the issue of the compatibility of ( 30b ) 



with (30a I. As has already been mentioned in section 
PH for p — > oo and \i = 0, both cases should match. A 
naive limit of (30b) shows that this is not obvious. The 



point here is that one sho uld c arefully look at the time 
dependence. Expressing (30b) in terms of the number 
of e-foldings, N := ln{a/ai}, shows indeed the desired 



coincidence. So for all plots related to power-law cases, 
we will make the replacement 



t = 



(32) 



A straightforward calculation shows the relation between 
t and N, 



N 



-E-[i+Mp}] =t + ln{p} + 0{p- 1 ). (33) 
p + 1 



This means that for large p, t = N up to a shift which 
originates from different time normalization of the power- 
law and the de Sitter case where t = N. This shift is 
indeed visible in figures [4] and [6j where only the last part 
of the transient phenomenon does show up, contrary to 
the corresponding exponential inflation plots, where a 
larger part can be observed (details below). 

In figures [l] (and U we show the effective long- 
wavelength power spectrum V as a function of k for fixed 
time Ht — 10 and mass /i = 0.1 (t = 4 and p — 12) with 
e = 10~ 2 and k — 10~ 3 for the de Sitter (massless power- 
law) model. One can see that it diverges stronger than 
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Figure 1: Power spectrum V of a massive test field (fi — 0.1H) 
for exponential inflation at tH = 10 as a function of comoving 



momentum k. The parameters of the filter function (27 1 are 
fixed to k = 1CT 3 and e = 1CT 2 . 
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Figure 2: Power spectrum V for a massless test field for power- 
law inflation (p = 12) at t = 4 [c.f. ( |32[ ) ] as a function of 
comoving momentum k (in units of l/£i). The parameters of 
the filter function ( 27 1 are fixed as in figure [T] 



the noiseless power spectrum as k tends to zero, putting 
therefore more correlation on large scales. One further 
sees that the part /c 3 ctG§ approximates the function the 
full power spectrum V in the infra-red as well as that 
the noiseless piece fc 3 Go gives a suitable ultra- violett ap- 
proximation. In all the plots we took the filter function 



(35a I 



Figure [3] shows the time behavior of the comoving k* 
for exponential inflation for different values of the mass /i. 
Figure [4] displays the same for the power-law model. The 



solid rays represent the analytic approximations (30a 
and ( |30b| , respectively, while the dashed curves arc the 



full results, obtained numerically from (25 1 using the pro- 
cedure FindRoot of Mathematica 6. Well below this 
borderline the two-point function obeys dimensional re- 
duction, while well above ordinary scaling holds. Af- 
ter an initial transient phenomenon, which' duration de- 
pends on the specific choice of e and k, the comoving 
transition scale decays exponentially fast. Hence, the 



Figure 3: Comoving transition scale k, [c.f. (30al] 
for exponential inflation as a function of cos- 
mic time t (in units of H) for mass /J./H = 
0.1, 0.2, 0.3, 0.4, 0.5, 0.6 (dashed lines, top to bottom). Dash- 
ed curves are numerical results, colored solid lines are 
analytic approximations, and enveloping black lines are | a(t) 
and the asymptotic f orm (31 1, respectively. The parameters 
of the filter function (27 1 are fixed as in figure [l] 
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Figure 4: Comoving transition scale k, for power-law in- 
flation as a function of modified cosmic time t [c.f. ( 32 1] 
for p = 5, 10, 12, 14, 16, 18 (dashed lines, top to bottom) . 
Dashed curves are numerical results, colored solid lines are 
analytic approximations. The parameters of the filter func- 
tion ( 27 1 are fixed as in figure [l] 



dimensional-reduction contribution is pushed to larger 
and larger scales as time increases. This therefore guar- 
antees that quantum modes induce only a minor change 
of the spectral index on sub-horizon as well as on mod- 
erate super-horizon scales. 

For concreteness, let us consider a mode with comov- 
ing k = 0.05 H. At time t = it is within the region of 
ordinary scaling, suffering at most slightly from dimen- 
sional reduction. This mode enters then, after roughly 
two e-foldings, the region of broken scale invariance, but 
leaves it at the latest (for /x = 0) after seven e-foldings 
and stays eternally in the scale-invariant regime, which 
itself grows exponentially fast. 
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C. Filter Functions 

The discussion of possible smooth filter functions and 
their influence on the phenomenon of dimensional reduc- 
tion shall be subject of this subsection. In particular, we 
study their effects on the transition scale fc*(t). 

Let be the Heaviside function. Then, we define a 
filter function as a function W K depending on a parameter 
k (controlling the width of the transition) such that 



lim W K = 6. 

K— >0 



(34) 



One may divide filter functions fulfilling (34) into two 
classes: Those for which — W K has an infinite support 
I and those for which I is compact. Let us now discuss 
those two cases separately. 



1. Non-Compact Support 
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Figure 5: Influence of filter function shape on the comoving 
transition scale fe* for exponential inflation as a function of 
cosmic time t (in units of H) for mass /J./H = 0.1. A smooth- 
ing k — 10~ 3 and short-wavelength cut e = 10 -2 are chosen. 
The variable k' is a short-hand notation for (fcr — e)/«. 



Some well-known smooth filter functions are: 



tan 



S(¥) 

1 

1 + e « 



1 

" 2' 
1 

2' 



1 
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tanh 
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(35a) 
(35b) 

(35c) 

(35d) 

(35e) 
(35f) 



Here, erf is the error function, defined by erf(z) := 
^= Jq dt e _ * and Si is the sine-integral function, defined 

by Si(z) := /*dt sm{t)/t. 

The following two figures show the dependence of the 
transition scale on various filter functions, where we 
choose the functions ( 35a I , ( 35d I and ( 35f I . It can clearly 



be seen that only the quantitative behavior changes, i.e., 
the position of the bump, which marks the end of the 
transient phenomenon, and not the qualitative shape. 
Note, that the curves are rescaled by a fixed factor (one 
for each filter function) and that those functions all have 
precisely the same asymptotic behavior. These state- 
ments concern both the exponential as well as the power- 
law case. 

In figure [7] we display the dependence of fc* on various 
values of the width parameter k. We find that decreas- 
ing this parameter shifts the curves downwards, there- 
fore pushing the dimensional reduction effect to larger 
and larger scales. This can also be seen directly from the 
late-time formulae (30a I. The power-law case behaves 
similar. 

An observation is that in the limit k — > the dimen- 
sional reduction part is absent. This is clear, because 
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Figure 6: Influence of filter function shape on the comoving 
transition scale fc* for power-law inflation as a function of 
modified cosmic time t [c.f. |32l] for p = 12. Filter argument 
and parameters are fixed as in Fig. [5] 
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Figure 7: Dependence of the comoving transition scale k„ 
on the width k of filter (271 for exponential inflation as a 
function of cosmic time t (in units of H) for mass j-i/H = 0.1 
and e = 10~ 2 . 
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there is no mixing of the short quantum modes with 
the long classical ones in the free (Gaussian) theory as a 
sharp cut off is introduced. In a slightly different setup, 
with a filter function having only one parameter, i.e., 
e = K, this has already been noted in reference [55] , 



2. Compact Support 



We note that all of the filter functions ( 35af f ) do not 
have a lower bound on their support — a crucial ingredient 
for the large occurrence of dimensional reduction in the 
far infra-red. This can be seen as follows: From (23) 
one observes that the dimensional reduction part, cGg, 
is proportional to <r, which itself is related to the filter 
function W K in such a way that if W K ^ 8 on an interval 
I only, cr = outside of I. Consequently, G = G , i.e., 
dimensional reduction is absent, on the complement on 
R\l. 

However, this does not mean that we can forget about 
dimensional reduction in the context of stochastic infla- 
tion: Any smooth filter function W K will definitely cause 
a deviation from scale invariance (see below) , although it 
might be that this deviation disappears for scales outside 
the support of W K . Further, the results of this work are 
not arbitrary, since an ultimate derivation of the stochas- 
tic inflation paradigm from first principles will single out 
a concrete filter function. 

We now study the effect of generic filter functions with 
compact support on the power spectrum. Since u is a 
free mode function, obeying (□ + H 2 )u = 0, it follows 
that the field equation for the long wavelength part does 
only contain derivatives on W K . For I = (—k,+k) one 
may choose 



(36) 



as a prototype filter function, where k! := (fcr — e)/K. 
Hence, its derivative W^ is given by 







k' 


< -K 


W K (fc') = < 


f 


k' 


e I 






k' 


> +K 







k' 


< —K 


W' K (k>) = < 




k' 


e 1 






k' 


> +K 



(37) 



leading to the approximation 

{0, k' < -K 

expll-^^}, fc'Gl . (38) 
0, k' > +k 

Figure [8] shows the influence of I ^ K on the power 
spectrum for the filter function due to ( 38 1 . Although 
for kr > e + k the dimensional reduction effect disap- 
pears, one clearly has an effect inside the interval I. As 
the size of I shrinks, the domain in wave-number space for 
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Figure 8: Effective power spectrum V(k) for filter ( |38[ ) with 
compact support I = (— k,-\-k) for different values of n. V(k) 
is evaluated for exponential expansion at t = 6.5H with /i = 
and e = 10~ 2 . 



which the dimensional reduction part is dominant does 
also size down, albeit the magnitude of the effect of wave- 
length separation increases considerably. This is reason- 
able, since the (step-function) limit k — > contains sec- 
ond derivatives on the Heaviside function corresponding 
to the pole forming for ( 38 1 . 

In figure [9] we display the effective power spectrum 
V{k) for different time. As anticipated, the dimensional 
reduction 'bumps' decline as time increases. One also ob- 
serves the same behavior as in figure [3] namely the grows 
of the comoving transition scale, up to the point where 
suddenly the dimensional reduction effect disappears (it 
is sub-domonant in the shaded region). Plot 10 visualizes 
the ratio 



n 



^Gj 
G n 



= aG n 



(39) 



of the dimensional reduction part to the noiseless part, 
evaluated at the most ultra-violet peak. One sees that the 
former diminishes exponentially fast in time. Hence, af- 
ter some few e-foldings, the transition region (from short 
to long wavelengths) becomes unimportant (shaded re- 
gion in figure |10| and the classical power spectrum pro- 
vides an excellent approximation. 



D. Modified Gaussian Fluctuations 

One may connect the replica structure a to a non- 
linearity parameter g NL , which shall now be defined via 

^(t, k) = <p?(t, k) - g nh (t, k)(<pf(t, k))\ (40) 

where (p G {t,k) is a free Gaussian field. On the level of 
propagators, this translates to 

G L (t,fc) = G L0 (i,£:) + 3 5NL (i,fc) 2 G L0 (i,fc) 2 (41) 

and hence 

a(t,k) =3<? NL (t,fc) 2 (42) 
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Figure 9: Effective power spectrum V(k) for various times 
and k = 10 -3 , otherwise as in Fig. [8] 
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Figure 10: Time dependence of the ratio rj of the dimensional 
reduction part to the noiseless part of the effective power 
spectrum for the filter ( 38 1 for exponential expansion with 
k=10~ 3 



6 = 10~ 2 and fi = 0. 



can be directly read off, using equation ( 23 1 . g NL mea- 
sures the influence of the quantum fluctuations, picked 
up by a smooth filter function. Formally, it resembles 
an effective non-Gaussianity parameter |34j for the long- 
wavelength modes. However, this association is mislead- 
ing since the theory we work with is Gaussian (but with 
non-trivial replica structure). 



1. Non- Compact Support 

Let us first consider the case of a non-compact filter 
function, which we choose to be (27 1. Figure 11 shows 
the dependence of g NL on the comoving momentum k for 
various values of p for fixed time t = 10H, using equation 
(42 1. Firstly, one sees that increasing p, shifts the curve 



upwards, and secondly, one observers a divergence in the 
infra-red displaying the effect of dimensional reduction. 
For k^$> k t one obtains a scale-invariant spectrum. Fig- 
ure 12 visualizes the same with p = for various p, where, 
t has been fixed to t — 4. One observes that increasing p 
lowers the curves with converge toward their asymptotic 
value for p — > oo, which is the same as the p — > limit 



Figure 11: Non-linearity parameter t/ NL for exponential in- 
flation as a function of comoving momentum k (in units of 
H) for mass p,/H = 0.1 (uppermost), 0.2, 0.3, 0.4, 0.5 and 
0.6 (lowermost) at Ht = 10. For the filter (27} we fix 
k = 10 -3 and e = 10 -2 . The crossings of the dotdashed black 
line with dashed lines indicates the value of fe» for different 
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Figure 12: Non-linearity parameter g NL for power-law infla- 
tion as a function of comoving momentum k (in units of 1/ii) 
for p = 10 (uppermost), 12, 14, 16, 18 (lowermost) and t — 4. 
Filter and parameters are as in figure [TT| 



of the exponential case as noted in section ITT] 



2. Compact Support 



To discuss the effect of a compact support I on the 
non-linearity parameter <? NL , we choose the filter corre- 



sponding to ( 38 1 . The dependence of g NL on the comov- 
ing momentum k is depicted in figure [l_3j All curves 
shown in this plot are strictly zero outside the plotted 
(momentum) intervall. This means that there is only a 
small (~ k) window around e in which the wavelengths- 
separation effects play a role at all. Albeit, for suffi- 
ciently early times, the transition effect becomes indeed 
pronounced. Figure [14] shows the time dependence of 
the non-linearity parameter <? NL for various masses. Af- 
ter some few e-foldings, g NL is completely negligible and 
Gaussianity of the fluctuations in the proper sense holds 
true. 
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Figure 13: Non-linearity parameter g NL as a function of co- 
moving momentum k (in units of H) for de Sitter infla- 
tion with mass Ll/H = 0.1 (uppermost), 0.2, 0.3, 0.4, 0.5 and 
0.6 lowermost) with Ht — 1. The filter ( |38| l with compact 
support is used with « = 10~ 3 and e = 10~ . 
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Figure 14: Non-linearity parameter g NL as a function of 
time t (in units of H) for de Sitter inflation with mass 
fi/H = 0.1, 0.2, 0.3, 0.4, 0.5, 0.6 (dashed lines, top to bottom) 
evaluated at the first ultra- violett bump of Fig. |13| Filter and 
parameters are as in Fig. ]13| 



their support. For filter functions with infinite sup- 
port we find the phenomenon of dimensional reduction 
on super-horizon scales. It heavily amplifies the power 
spectrum in the infra-red. The time-evolution of the 
long-wavelength field pushes the dimensionally reduced 
region exponentially fast to unobservable scales. Taking 
the limit of vanishing width of the filter function n, i.e., of 
a sharp separation of long- and short-wavelength modes, 
has the same effect. 

For filter functions which deviate from the step func- 
tion on a compact intervall I only, we show that the 
smooth separation might also lead to strong modifica- 
tions of the power spectrum. However, this effect is lim- 
ited to I and decreases exponentially fast in time, becom- 
ing negligible after a few e- foldings. 

Our findings provide further support for the self- 
consistency of the idea of inflation. Either regions of 
broken scale-invariance with extraordinarily large fluc- 
tuations disappear faster than any causal patch of the 
universe expands (non-compact I), or large extra power 
is strongly damped by the time evolution (compact I). 

The huge effect due to quantum noise on large super- 
horizon scales (may they occur on a finite or infinite mo- 
mentum range) do not permit us to speak about the spec- 
trum of fluctuations in the usual, perturbative sense and 
clearly signal a breakdown of ordinary perturbation the- 
ory. This noise modifications further display the failure 
of the test-field assumption, since in the situation at hand 
it is no longer valid to neglect the back-reaction of the 
field on the geometry. In the case of compact I it is pos- 
sible to avoid the breakdown of perturbation theory and 
the test field assumption by an appropriate choice of the 
filter width k. 

The next step will be to apply our new methods to 
interacting salar fields and, in turn, discuss replica sym- 
metry breaking solutions. Another open task is to study 
more general space-times, such as those modelling inho- 
mogeneities. 



V. SUMMARY AND CONCLUSION 

In this work we have studied the large-scale behav- 
ior of the power spectrum of the long-wavelength part 
of TV-component scalar test fields in curved space-time, 
using a stochastic description for the quantum modes, 
introduced by Starobinsky pQ. We focus on the two im- 
portant cases of a spatially-flat Friedmann geometry with 
an exponential and a power-law scale factor. The effec- 
tive spectral index n v is calculated in the framework of 
replica field theory a la Mezard and Parisi [23], which we 
recently introduced in a cosmological context [22]. Using 
a Gaussian variational approximation |35j we derived an 
expression for the physical propagator G and thus for the 
power spectrum of long- wave length fluctuations. These 
methods allow us to study the spatial behavior of arbi- 
trary long- wavelength two-point correlation functions. 

A discussion on possible filter functions has been given 
with special focus on the aspect of the compactness of 
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Appendix A: NOISE DISTRIBUTION AND 
REPLICA TRICK 

In this appendix we present some mathematical details 
and methods needed for our analysis. Because those tech- 
niques originate from the study of disordered systems in 
the framework of statistical field theory, we Wickrotate to 
an Euclidean action. All formulae can trivially be trans- 
formed to those for spaces with signature (— , +,...,+) 
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by simply performing a Wick rotation in the temporal 
component. Further we generalize from the free case, 
discussed in the main text, to general stochastic interac- 
tions [c.f. ( 20 1] and work in d space-time dimensions. 

First, we explain the replica trick, which provides a 
simple technique to easily compute stochastic averages 
of the generating functional, i.e., general n-point corre- 
lations. Then, we study the noise distribution of the 
involved random variables, focussing on a Gaussian one. 
This has the advantage that one can obtain results ana- 
lytically, being exact in the case of a free theory, though 
in general only approximative. 



1. Replica Trick 

Observables in stochastic systems, i.e., those carrying 
random variables with a certain probability distribution, 
adhere additionally to the classical average over (p L — 
an average over the noise variables h, which shall here 
be denoted by a bar. In particular, to compute corre- 
lation functions and hence power spectra, one needs to 
compute the noise average of the generating functional 
of connected n-point functions (in the context of statisti- 
cal physics, the object of interest is the free energy) and 
thus ln{Z[f]}, for an external current ]. This is difficult 



to achieve directly since one has to average a logarithm 
of a path integral over an exponential. Instead, one uses 
the so-called replica trick |23j : 



5" 



5jh On 
= lim 



m->o m8] ll (x ll ) ... Sj in (x in ) 

- lim 1 — 

m^o m 5j n (x n ) ... Sj in (xi n ) 

- lim 1 &n 

m->o m Sj i± (x i± ) . . . 6j in (x in ) 



Mm} 

ln{l + mln{Z[J]}}} 



ln{exp{m ln{Z[f]}}} 

, ln(z™[J]} . (Al) 

m^o m dji^Xij . . . oj in {x in ) I J 



lim 



6" 



This means, that one first has to compute Z m for m 
integer. Then, if the result is an analytic function of 
to, one performs an analytic continuation and takes the 
limit to — > 0. Z m , describes to non-interacting random 
systems (replicas). By performing the average Z m , this 
to interacting systems become coupled by means of the 
stochastic average. 

Let us now present, as an application of the replica 
trick, the computation of the full two-point function. 



G(x,y) 1= - (0{x) ■ 0(y)) ~ - ^- ^- 



m{Z[j]} 



7= 



iAJl hi * hW 

m^o Nm Sj(x) ■ 5j{y) I L J J 



lim 

m— >o Nm 



fWi I II v ^p^ Y ■ Vd(y) expj - ^2s[0 b ] > 

J a=l c,d=l { 6=1 J 



=2/11 V ^a} Y Vc{x) exp< -^2 S[<p b ] •/[] V[0 b ] <Pd(y) expl - ^ S[<Pb] \ 

* m •> a=l c=l I 6=1 ) J 6=1 d=l I 6=1 J 



lim 



1 1 



+0 Nm Z r 



■ / n y ^ x "> ■ exp j ~ y s ^ 

J o=l c=l ^ 6=1 



1 1 r 

lim — y^G aa (x,y)= lim — Tr (G ab (x,y)) 

m^O m £ — • m^O TO L 



(A2) 



with Z := Z [ j* = ] and 



G ab (x,y) := 



N 



1 /* ( 

= lY[V[<Pc]0a(x) ■ b (y)expl - ^2S[0 d ] 

* J c=l L d=l 



(A3) 



In the second step of the first line in ( A2 1 we account for the infra-red parts only, where we assume that 



(<p(x)) ■ (<p(y)) 



(A4a) 
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has a milder infra-red behavior than 



(0(x) ■ 0(y)) 



(A4b) 



i.e., we demand that for k — > the disconnected and the 
connected correlation functions coincide, 



(<p(x)) ■ (<p(y)) 
(0(x) ■ <p(y)) 



(A5) 



This is exact for the free case discussed in the main text 
as well as for an even action, i.e., S[<p] = S[— <p]. 

2. Noise Distribution 

Here, we present some general considerations on the 
distribution for the random variables h^ and introduce 
some of the statistical physics jargon. For simplicity, we 
shall assume zero mean, i.e., /i\ = 0. This does not 
restrict the following discussion, and since we will only 
average exponentials, we can always shift fi to zero with 
the effect of adding a potential, which' coefficients are 
just given by the corresponding average values, i.e., h\ <-> 
H\. We will see later on that such potentials generate 
mass corrections. The second cumulant of V shall then 
be taken of the form jJD] 

^a b{x,y) 
N 

'(A6) 



V{$ a {x),x)Y{0 b {y),y) = ^(x,y)NR 



where, for later convenience, we rescale by the number 
N of field components. The function R reflects the cor- 
relation in replica field space, with, e.g., 



S a b(x,y) : = 0a{x) ■ tp b (y) 



Zab(x,y) ■= [0a(x) ~ (Pb{y)\ 



(A7a) 



(A7b) 



We refer to (A7a) as product correlation, while (A7b) 



is called difference correlation. In the following we will 
mainly use the former (c.f. the discussion at the end of 
this appendix). Of course, the concrete form of the func- 
tion R, as well as of its argument have to be determined 
from first principles. The space-time correlation cj>(x, y) 
is called short-range, if <j>{x,y) = p(t,t')5 d ~ 1 (x — y), and 
long-range for all other cases (c.f. appendix |B 2| and also 
[36]). Note, that for our example of a free field, the Gaus- 
sian distribution is exact and that the function R is lin- 



ear. The difference between (A7al and (A7bl is that the 



latter has the so-called statistical tilt symmetry: 

<fi a (x) -> a (x)+g(x), 



(A8) 



where g(x) is some function without (replica) index. 

When performing the noise average over correlation 
functions, one has to average the measure, and thus, 



after replication, exp{— ^ Q Assuming product- 
correlation and the matrix (Aa,a') to be diagonal in com- 
ponents and position (short-range correlation), the part 
containing noise may be calculated as 



exp< 



h A (fx-.. <p\ 



cx V 



[{h}] cxp<j ^h X ifix ■ ■ • ipx - i ^h A A A .A'h A ' 



A. A' 



exp< ^^(^a-.-^aA x x,\'tp\ ■ ■ ■ Vx> \ 



(A9) 



0a(x) ■ Ifb(x)' 

N 



where J := J d d x. 

This is perhaps a good place to mention a difference 
between what is described here and what is described in 
the field of disordered systems. In the latter one studies 
macroscopic objects, e.g., a crystal, with defects, e.g., 
vacancies, wrong atoms or molecules, misaligned layers or 
substrate impurities. These kinds of disorder is mimicked 
by some random variables - a prime example of such a 
system is the random field Ising model |37j . 

Averages over disorder are thought of as averages over 
different realization, i.e., practically different pieces of a 
crystal, for instance. Mezard and Parisi [33] made a con- 
sistent replica field theoretic approach to those systems, 
which' m copies, arising due to the application of the 
replica trick [c.f. (All], are viewed as respectively differ- 
ent, coming from different realizations, carrying different 
random variables, with distinct correlations among them. 

This is, however, not the case in stochastic inflation, 
where the noise arises from the short wavelengths of some 
quantum field within one and the same system. For the 
case under consideration the replica trick (All gives us 



a simple procedure to perform stochastic averages over 
the generating functional. There, each Z[f\ is a function 
of the random variables {h}, so its averaged m'th power, 
Z m [j], is too. It is only the integration variables cp L that 
acquires an additional label, namely the replica index 
a = 1, . . . , m. 



Appendix B: VARIATIONAL METHOD 

This appendix is devoted to the application of the 
Feynman- Jensen inequality and a Gaussian variational 
principle to derive variational equations for a variational 
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propagator G ab , which allows us to obtain a closed ex- 
pression for the power spectrum of the long-wavelength 
modes. 



1. Feynman- Jensen Inequality 

To perform the stochastic average over the (quantum) 
noise variables {h}, we use the replica trick (c.f. appendix 
A 1 1 and obtain 



/m ( rn 

n^lexpi -J2S[0b,j\ 
a=l I b=l 

rn 

Y[V[<p a ]exp{-S( m \<pj]}, 



(Bl) 



with the replicated action (for short-range product- 
correlation, for simplicity) 



-k) 



Z a,b=l Jx >V 



N 



(B2) 



Here and in the following, we refrain from writing the 
different replica fields in the argument of the action, i.e., 



[0] ■= [01, ■■■,07, 



(B3) 



for simplicity. 

We proceed with a Feynman- J ens en variation princi- 
ple |35j and therefor define the Gaussian trail action 



So[0] = 2 E 



,6=1 



t.t'k 



G- l ab {t,t',k)0 a {t,k)-0 b {t',-k), 

(B4) 

where we make the following ansatz for the propagator 
G _1 ab := (Gq 1 + a c + a aa ) S ab - a ab . (B5) 



Let us comment on the structure of ( B5 1 . On the main 



diagonal (in replica space) we have the inverse of the 
noiseless propagator Gq 1 plus some mass correction cr c , 
to be determined later, e.g., by the variational principle 
described below in this section. This alone would not 
only be trivial but also inconsistent, as we will see later. 
Hence, the off-diagonal part is filled by some, a priori 
unknown, replica structure o~ abl which, in general, can be 
time dependent, and, if one includes long-range noise cor- 
relation (c.f. appendix B2|, also momentum dependent, 
directly affecting the scaling behavior of the power spec- 
trum. Thus, although this variational method only gen- 
erates a self-energy contribution, its off-diagonal replica 
structure might have a viable influence on large-scale cor- 
relations. 



The Gaussian variational method becomes exact in the 
limit N — > oo and allows one to go beyond ordinary per- 
turbation theory. It is based on the following Feynman- 
Jensen inequality [35) 

ln{Z} > ln{Z } + (s^ m) - S(" l >) Q , (B6) 

where the subscript refers to the variational action 



(B2) and we temporarily Wickrotate to Euclidean sig- 



nature. Equation (B6l can easily be proven by using 
the Jensen inequality exp{(. . .)} < (exp{. . .}) [3H], which 
comes from the convexity of the exponential. The prob- 
lem is to find the best G La b, i.e., the best a ab , satisfying 
(B6) by maximizing the r.h.s. of (B6l. 



Computing T ml per component and spatial volume 
yields 



1 m 



t.t'.fc 



G - 1 (t,t',fc)G QQ (t,i',fc) 



iVVol(d- 1) 

-\\ Trln{G(M',fc)} + C (B7) 

2 Jt,t'.k 



,6=1 



cj>(t,t')R[ / G ab (t,t',k) ) , 

t,t' \Jk 



where we temporarily switched to finite spatial volume 
Vol(e? — 1). The constant C, which vanishes after varia- 
tion, includes (So) as well as terms from J-q. 

For difference-correlation (A7bl one just has to replace 
the argument of R to 



(G„„(t, k) + G bb (t, t',k)-2 G ab (t, t', k)) . (B8) 



The "hat" over the function R in (B8) is defined through 

R((.)o) :=<R(-))o- (B9) 

In the limit N oo, averaging and applying the (ana- 
lytic) function R commute and so we drop the hat when 
such a limit is considered. 



The variation of (given in (B8l) with respect to 
the m 2 variational parameters G a b gives for a ^ b 



<Tab{t) =<K*)R' / G ab (t,k) 



and 



<r c (i) = -<£(*) R / G aa (t,k) 



where we use 

SGgb(t,k) 

SG cd (r,p) 



= 5(t-r)5 {d - 1 \k-p) 6 ac 5, 



(BlOa) 



(BlOb) 



(Bll) 



and define the equal-time entities s ab (t) := s ab (t,t), 
s c (t) := s c (t,t), G ab (t,k) := G ab (t,t,k), and G (t,k) := 
G (t,t,k). 
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Again, for difference-correlation we obtain a similar 
result as (BlObl (c.f. equation (3.12) in [23]) with the 



argument of R' replaced by 



(G aa (t, k) + G bb (t, k) - 2 G ab (t, k)) (B12a) 



plus a global minus sign. Additionally, we find in this 
case 



a c (t) = -}^<J ab {t). 



(B12b) 



a=l 



It is important to note that any interaction not con- 
taining noise just modifies the self energy through its di- 
agonal structure in replica space and that the free Gaus- 
sian case studied in the main text, necessarily yields 
replica symmetry. 

The physical interpretation of the saddle-points equa- 
tions |B5]), flB10a| b) and ( |B12a| b) is the following: The 
replica structure a is a generalized self energy, i.e., the 
sum of generalized tadpole diagrams (c.f. the discussion 
in section 3 of [23] ) . This follows from an expansion of the 
stationarity equations ( B10a[ b) and ( B12a[ b) in powers 
of(G„ 6 ). 



2. Long-Range Correlation 

We now consider the case where the noise-correlation 
is non-local, or of so-called long-range-typ e [55] , This is 
described by the following correlation c.f. ( A7a b)] 



V{<p a (x),x)V(<p b (y), y) = <f>(t,t / ,x-y)Nx 

'Z ab (x,y)\ (B13) 



x R 



N 



Going through analogo us step s as in the previous section, 
we obtain pendants to ( B10a[ b) (product correlation), 



<r ab (t,p)= / 4(t,x)e- i r x R'[ / e- lkx G ab (t,k) 



(B14a) 



a c (t, P) = -J 4>{t, x) e-*- R'U e- ik x G aa (t, k) 

+ 2U'(jG aa {t,k)), (B14b) 

and to ( B12a| b) (difference correlation), 

a ab (t,p) = -2 [ <P(t,x)e- [ P x R'f f C - ikx (G aa (t, k) 

J x \J k 

+ G bb (t,k)-2G ab (t,k)Yj , (B14c) 



°c(t,p) 



E 

«=i 

a#6 



Vah{t,p), 



(B14d) 



where a ^ b for (B14al and (B14cl. We introduce a non- 
random potential U (with suitable rescaling by factors 
of N), which is so far taken to be arbitrary. It may 
arise from non-zero averages of the random distributions 
(c.f. the discussion at the beginning of appendix A 2 1. 
The last equation is restricted to U' = 0. 



Equations (B14a d) already contain our main result, 
namely the dimensional reduction on large scales. It also 
shows that the replica matrix is in general space-time 
dependent, which affects the scaling behavior of the two- 
point function (c.f. 32J). Clearly, the short-range varia- 
tional equations are included in (B14a d). 



3. Replica Symmetric Propagator 

In the previous section we derived the variational equa- 
tions (B14a d) for a general matrix (a ab )- It is impor- 



tant to try the simplest ansatz, which consists of taking 
a ab = a for all a ^ b, meaning that different replicas cou- 
ple all in the same way among each other. This replica 
symmetric case is exact for the free case discussed in the 
main text of this article. One finds 



(G^)- 1 - (Go 1 + <r c )l - a(3 - 1) 
- (Gq 1 + (7 C + cr)l - crJJ 



(B15) 



Gn 1 ! -o-JJ. 



Thus, the inverse G ab (t,k) has the form 

(G ab )(t, k) = G (i, k)l + a(t, k) Go(t, fc) 2 J, (B16) 



where the physical limit m — *■ has been taken and we 
define the m x m-matrix JJ by JJ.y = 1 for all It 
has the property Jl 2 = rraJJ, obtained by inverting ( B15 1 , 



is equal to zero in this limit. We observe that the limit 
of vanishing correlation, i.e., a — > 0, gives back the free 
power spectrum as expected and obtained in [33 , with 
rather different methods. 

The physical propagator G(t, k) of the long- wav elen gth 
field is obtained from G ab (t,k) via [c.f. equation (A2)] 



G(t,k)= lim -Tr[(G Qb )(i,fc)], 



(B17) 



which is simply the arithmetic mean of the trace of the 
replica matrix propagator. 
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